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This paper is concerned with utilizing the acoustic analogy approach to predict 
the sound from unheated supersonic jets. Previous attempts have been 
unsuccessful at making such predictions over the Mach number range of 
practical interest. The present paper, therefore, focuses on implementing the 
necessary refinements needed to accomplish this objective. The important 
effects influencing peak supersonic noise turn out to be source convection, mean 
flow refraction, mean flow amplification, and source non-compactness. It appears 
that the last two effects have not been adequately dealt with in the literature. The 
first of these because the usual parallel flow models produce most of the 
amplification in the so called critical layer where the solution becomes singular 
and, therefore, causes the predicted sound field to become infinite as well. We 
deal with this by introducing a new weakly non parallel flow analysis that 
eliminates the critical layer singularity. This has a strong effect on the shape of 
the peak noise spectrum. The last effect places severe demands on the source 
models at the higher Mach numbers because the retarded time variations 
significantly increase the sensitivity of the radiated sound to the source structure 
in this case. A highly refined (non-separable) source model is, therefore, 
introduced in this paper. 


1. Introduction 

The impracticality of directly calculating the sound field from the full Navier- 
Stokes equations creates a need for a viable reduced order model for these 
equations-typically referred to as an acoustic analogy following an approach 
initiated by Lighthill (1952). This is analogous to the situation in turbulence 
modeling, which is almost always based on some form of the filtered Navier- 
Stokes equations such as the Reynolds averaged Navier-Stokes (RANS) 
equations. While there is some disagreement about the proper choice of 
turbulence modeling equations, the situation in Aeroacoustics is even more 
contentious. Here, there is considerable disagreement about the appropriate 
starting equations. Perhaps because there is no single set of equations that is 
optimal in all situations. There does, however, seem to be a consensus about 
some of the requirements for such equations. First, they should be derivable from 
the Navier-Stokes equations and second, they should be formally linear (Dowling 
& Ffowcs Williams, 1983, p.157). 
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It, therefore, makes sense to describe this phenomenon by dividing the flow 
variables into their mean and fluctuating components. The Navier-Stokes 
equations can then be rewritten as a set of mean flow equations plus a formally 
linear set of equations for the fluctuating component of the motion. And, it turns 
out that the latter can be put into the form of the linearized Navier-Stokes 
equations (about the mean flow) but with source terms whose strengths can be 
represented by a four dimensional stress tensor, whose first 3 dimensions 
correspond to the usual fluctuating Reynolds stress with the additional dimension 
being associated with the stagnation enthalpy fluctuations. The true non-linearity 
of these equations is hidden in the source terms as well as in the non-linearity of 
the dependent variables. The latter causes no particular difficulty in the present 
context because one of these variables (the pressure-like variable) reduces to 
the ordinary pressure fluctuation in the far field. 

The apparent linearity can be exploited by using a Green’s function 
approach to separate out the so called “propagation effects” from the unsteady 
source fluctuations, which are represented by a generalized fluctuating stress 
tensor in the present approach. The result can then be used to express the far 
field pressure autocovariance as the convolution product of a “propagator” 

(which can be calculated from the Green’s function) with the two point time- 
delayed correlation of the fluctuating stress tensor-which is a statistical entity of 
the type that is typically measured in turbulent flows. It is also the quantity that 
ultimately has to be modeled in the present acoustic analogy approach. It is 
expected that most of the non-local “propagation effects”, which would be very 
difficult to distinguish from the turbulent fluctuations-and, therefore, very difficult 
to model- have been removed from these stresses. An attractive feature of the 
original Lighthill (1952) approach is that the source strength is represented by a 
single stress tensor, which insures that the far field pressure autocovariance can 
be expressed in terms of a two point time delayed correlation tensor. This 
significantly simplifies the interpretation of the source (Ffowcs Williams, 1963, 
1969) and thereby aids in the construction of appropriate models for this quantity. 
The present result (which is similar to the one in Goldstein (2003)) also has this 
important property which, together with the two requirements set out at the 
beginning of this section, effectively restricts the form of the analogy to relatively 
trivial rearrangements of a single set of equations. 

Current state of the art noise prediction methodologies, typically, use 
empirical models for the unknown Reynolds stress correlations and RANS based 
approaches (typically of the k-epsilon type) to calculate the mean flow and mean 
turbulent kinetic energy. The former is then used to calculate the coefficients in 
the Green’s function equations, while the latter is used to determine the 
parameters empirical models for the fluctuating Reynolds stresses. 

This approach is very general and should, in principle, apply to any 
turbulent flow. But jet flows are nearly parallel, and the mean flow is usually 
approximated by a much simpler unidirectional transversely sheared flow. The 
resulting acoustic equations can then be reduced to a single 3 rd order equation - 
frequently referred to as Lilley’s (Lilley, 1972, 1974) equation. 
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This type of model appears to be adequate for subsonic, unheated jets 
(Khavaran & Bridges, 2004) but becomes unphysical at supersonic acoustic 
Mach numbers where it predicts infinite acoustic pressures. The difficulty is due 
to the breakdown of the parallel flow model and can be eliminated by retaining 
the correct non-parallel mean flow. But this would greatly increase the 
complexity of the Green’s function computation. 

The present paper introduces a perturbation approach that takes 
advantage of the small jet spread rate. It uses matched asymptotic expansions 
to obtain a locally parallel “outer solution” that applies in the main part of the jet 
where the parallel flow “propagator” is non-singular and a non-singular “inner 
solution” that accounts for the nonparallel flow effects in the vicinity of the so 
called critical layer where the outer (parallel flow) solution breaks down. The 
results are then combined to obtain a uniformly valid “composite” solution that 
produces physically reasonable results and is still not much more complicated 
than the parallel flow solution. The resulting “propagator” still becomes large 
when the source point is at the critical layer but, unlike the parallel flow result, 
remains finite there. 

A similar weakly non-parallel flow approach was used by Goldstein and Leib 
(2005) to construct a causal Green’s function that remains finite in the small 
spread rate limit. It results in an additional contribution to the usual parallel flow 
Green’s function that involves the linear instability modes of the jet. It was 
incorrectly stated in that paper that the result, which did not account for the 
critical layer singularity, was uniformly valid everywhere in the flow. The complete 
uniformly valid causal solution can only be obtained by combining the solution in 
the present paper with the Goldstein and Leib (2005) result. But this would be 
extremely complicated and our most recent computations show that the instability 
wave contribution tends to be very small at the relatively low acoustic Mach 
numbers being considered in the paper. We, therefore, decided to neglect that 
contribution (in the manner discussed in section 5) in order to focus on the 
critical-layer effects without introducing undue complication. 

While strong streamwise and relatively weak transverse coherence effects 
preclude the possibility of assuming that the sources are completely compact 
(i.e., of neglecting all variations in retarded time across the sources), it is still 
highly desirable to model the situation as closely as possible. The relevant 
source model should, therefore, account for streamwise variations in retarded 
time and be general enough to account for the long range turbulence correlations 
at the end of the potential core which, in particular, precludes using such 
assumptions as local isotropy and quasi-normality. The model should also be 
simple enough to produce easily computed expressions for the spectral tensor 
components. The present paper introduces such a model. Its most general form 
involves a large number of adjustable constants, which can be parameterized 
and subsequently determined from an auxiliary flow calculation. Current state of 
the art noise prediction codes, such as JeNo (Khavaran & Bridges, 2004), usually 
use steady k -s based computations for this purpose. But codes of this type only 
provide enough information to determine a small number of these parameters. It 
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is, however, highly likely that higher fidelity methods such as full Reynolds Stress 
Modeling or even hybrid RANS / Large Eddy Simulation approaches will replace 
the steady k -e computations. These codes can provide enough information to 
determine many more of the adjustable constants. 

Previous attempts -typically based on much simpler flow models — have not 
been successful at predicting jet noise over the parameter range of practical 
interest. They are usually unable to adequately predict the sound radiation from 
heated jets or the peak radiation from all supersonic jets. This has even caused 
some writers (see for example, Viswanathan, 2007) to recommend that the 
acoustic analogy approach be abandoned in favor of empirical correlations of the 
jet noise data base. But this approach is unlikely to be very robust beyond the 
parameter range of the available data. The present paper, therefore, focuses on 
implementing the necessary refinements into the acoustic analogy. It attempts to 
demonstrate that good predictions can be obtained when a high level of rigor is 
maintained by minimizing the approximations and keeping the formulas exact 
until the last stages of the analysis in order to make the required approximations 
as consistent as possible. We have chosen to not consider heated jets because 
the available computer codes for producing required information about the flow 
have not been adequately calibrated with an appropriate experimental data base. 
The focus is, therefore, on unheated supersonic jets, but the general formulas 
are applicable to heated jets as well. 

The important effects influencing peak supersonic noise turn out to be 
source convection, mean flow refraction, mean flow amplification, and source 
non-compactness. It is our contention that the last two effects have not been 
adequately dealt with in the literature. The first of these because, as originally 
noted by Phillips (1960), the usual parallel flow (i.e. Lilley equation) models 
produce most of the amplification in the so called critical layer where the solution 
becomes singular and causes the predicted sound field to become infinite. We 
deal with this by introducing a new weakly non parallel flow analysis that 
eliminates the critical layer singularity. This has a strong effect on the shape of 
the peak noise spectrum. The inability of previous attempts to deal with the last 
effect can be traced to inadequate source models. Retarded time variations can 
significantly increase the sensitivity of the radiated sound to the detailed source 
structure. A much more refined (non-separable) source model is, therefore, 
developed in section 6.2 of this paper. 

The overall plan of the paper is as follows: First, the basic acoustic analogy 
equations are set out in section 2 and a formal Green’s function solution is 
written down in section 3. The resulting expression for the pressure-like variable 
is an integral of a tensor product of a “propagator” that depends only on the 
mean flow variables, with a generalized four dimensional stress tensor that 
completely characterizes the turbulent velocity/enthalpy fluctuations. The 
propagator depends on the Green’s function for an arbitrary mean flow. The latter 
is simplified in section 4 by assuming that the mean flow is weakly non-parallel. It 
is shown that the strictly parallel flow result breaks down in the so called critical 
layer and matched asymptotic expansions have to be used to obtain a uniformly 
valid (i.e. non-singular) result. 
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The Green’s function solution of section 3 is used in section 5 to obtain an 
expression for the far field pressure autocovariance in terms of a two-point time 
delayed correlation of the turbulent velocity/ enthalpy fluctuations. This exact 
result, which shows that the far field pressure autocovariance can be expressed 
as a convolution product of a new propagator with this generalized correlation 
tensor, is the fundamental starting equation of this paper. It is used to obtain an 
expression for the far field acoustic spectrum in terms of the weakly non-parallel 
flow Green’s function constructed in section 4, which is then simplified by 
assuming that the transverse correlation length of the turbulence was small 
compared to the transverse length scale of the jet. The resulting expression for 
the acoustic spectrum is now purely algebraic and involves a product of a 
propagator and a generalized spectral tensor that accounts for the turbulent 
velocity/enthalpy fluctuations. It is only at this point that the results are restricted 
to “cold” jets by neglecting the enthalpy fluctuations. 

The propagation tensor is completely characterized by the mean flow which, 
with the present approach, is to be calculated from a RANS solution. But the 
generalized spectral tensor requires detailed information about the turbulence 
statistics, which can not be directly obtained from the steady RANS solution. It, 
therefore, has to be modeled, which is done in section 6. First “kinematic models” 
which exploit expected symmetry properties of the turbulence correlation 
functions, are introduced in section 6.1 to reduce its large number of independent 
components (45 in all, even if the enthalpy fluctuations are neglected). Then a 
very general correlation function model is introduced in section 6.2 to represent 
the spectra of the surviving components of this tensor. It involves a large number 
of parameters that can be related to the output of the RANS solution by making 
appropriate modeling assumptions. But before doing this, the general results are 
specialized to a round jet in section 7. The modeling assumptions are introduced 
in section 8 and the results are then compared with detailed jet noise 
measurements taken over a number of years at the Glenn Research Center. 

2. Basic Equations 

Assume that the pressure p , density p , enthalpy h , speed of sound c and 
absolute temperature T satisfy the ideal gas law equation of state 

p = pRT, h = c p T = c v c 2 / R (2.1) 

with R = c p -c v being the gas constant, and c p and c v the specific heats at 

constant pressure and volume, respectively. Then the Navier-Stokes equations 
can be written as 

Dp = 0, (2.2) 

D ^ +d irl^ <2 ' 3) 
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' ?p + p v ^ 


dp d 


\ y - 1 2 ) dt dx\y - 1 


cr„ 


■ + o-,.v. 


(2.4) 


where i,j = 1 , 2 , 3 , the operator D is defined by 


D f = $L+JL( v .f) 
dt dx ; v J J ) 


(2.5) 


for any function /, t denotes the time, jc = {x 1 ,x 2 ,x 3 ) the Cartesian coordinates, 
v = {v 1 ,v 2 ,v 3 } the fluid velocity. cr,j denotes the viscous stress, a 4i = -(y-1)#,-, 
the scaled heat flux vector, and y = c t c v the specific heat ratio. 


By introducing the new pressure like variable 


Pe=P + 1 ^ L ?V 2 


( 2 . 6 ) 


the Favre (1969) averaged Navier-Stokes (RANS) equations can be written in 
the virtually identical form 


D oP = °> 


(2.7) 


-tp 


( 2 . 8 ) 


D 


f YPe + P_~ 2^ ^ 


y - 1 2 


J 


dp^ _ _d_ 

dt dx 


f V - ' 


y-l 


j 


(2.9) 


where 




( 2 . 10 ) 


for any function / , the over bars denote the usual time average 


1 T 

s 1 •(*>')<*> 


( 2 . 11 ) 


with the dot being a place holder for p,v,,p, h, etc. and 
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• = (p*/p) 


( 2 . 12 ) 


denotes a Favre (1969) average, it being understood that the time derivatives 
drop out of the RANS equations (2.7) to (2.9), 


e U = -p K v 'i + 5 v ; 


fc-^pv' 2 +[a v; +(y-l)6 v4 a //t v;] 


(2.13) 


with v=l,2,3,4, 


v 4 = (r ~ 0 


h' + -v' 2 

v 2 , 


2 V , ( y -!).,2 


= + 


(2.14) 


and denoting the Kronecker delta in the usual notation. 
The primes denote the fluctuating variables 


p' = p-p, p' = p-p, h' = h-h, v'.=v i -v i , 


(2.15) 


etc., which upon introducing the (non-linear) dependent variables 


Pe=P' + ^2^\ PV' 2 - P V' 2 


(2.16) 


u i =A v i> 


(2.17) 


can be shown to satisfy the five formally linear equations (Goldstein, 2003) 


D °p' + — Uj =0, 

^ , dv,. d , p' d ~ d „ 

D u +u 1 p u = e 

0 ' J dXj dx, e p dXj 9 dXj lJ 


and 


(2.18) 

(2.19) 


^V c2 “ y+ ( rl ) 


, », 8e ij 


V 


dx. p dx 


j J 


d , „ dv, 

= —e tJ+ ( r pe,—, 


( 2 . 20 ) 


where 
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( 2 . 21 ) 




is the total mean flow stress tensor, 

? = ri>lp ( 2 . 22 ) 

the square of the mean-flow sound speed, and the source strengths 
e^j, v = l, 2, 3, 4 , which are given by 




= P . — P . 


(2.23) 


have zero time average. 

This shows that the entire acoustic source strength is given by the 
generalized (fluctuating velocity — sound speed) stress tensor e^j, v = i,2,3,4 ~ 

which is an important advantage of the present formulation because it puts all of 
the modeling on the same basis and thereby reduces the probability of 
over/under estimating one component of the source relative to another due to 
differences in modeling requirements (see paragraph following equation(5.29)). 

Notice that (2.7) to (2.9) will constitute a closed system of five equations in 
the five unknowns p,v.,^ that can be solved independently of the relation (2.6) 

between the effective pressure p e and the thermodynamic pressure ponce a 
particular turbulence model has been introduced. This is also true for the linear 
system (2.18) to (2.20) in the variables p',u i ,p' e — with their nonlinear relation to 

the physical variables being largely irrelevant because the variable p' e 

reduces to the acoustic pressure fluctuation p' in the far field where the sound 
field is to be calculated. 

3. Formal Green’s Function Solution for the Far Field Pressure 

The particular solution to equations (2.18) to (2.20) can, as noted in the 
introduction, be expressed in terms of the vector Green’s function (Morse and 
Feshbach 1953, pp. 878-886) g V(y (x,t\y,t) , which satisfies 


dv j dg 


D o giv+gj<,-^ + 


4(7 


1 d0„ 


dx i p dx. 


g5a= S i'j S ( X ~y) S ( t - Z ) 


(3.1) 


D og4a+(r- 1 ) 


dv t 

dx, 


S4c 


p dx f 


gu 


dc 2 g: 


JO- 


dXj 


= S, <T S(x-y)S(t-z) (3.2) 
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( 3 . 3 ) 


D ° 85a + dxj Sja °’ 

for a = 1,2,..., 5 and can be inserted into Green’s formula to obtain the following 
expression for the pressure like variable p' e (Goldstein, 2006) 


p'e( x ’ t ) = -\\ y ^( x Ay^) e lj{y^) d y dz ’ ( 3 - 4 ) 

V —00 

where V denotes integration over all space ,the propagator y k/ is 

defined by 


/ j ^ d S^(x,t\y,T) dv 

M y>*) s ( y )^r g44 M* T )* 


(3.5) 


v 4 = 0 , and we have integrated by parts to transfer the derivatives from the 
source term tog 4v . 


4. Weakly Non-Parallel Flow Approximation 

As noted in the Introduction, the propagator y (j,T|jc,r)can be greatly 

simplified by assuming that the mean flow is locally parallel. In fact, Appendix A, 
which describes the mean flow perturbation expansion, shows that equations 
(3.1 ) to (3.3) reduce to the inhomogeneous Rayleigh equations (i.e., the 
equations governing the inviscid instability of a parallel flow) in this limit 


%i 8ia + g4a + 5/1 ^ g j° = 5ia5 (* ' “ J ) ' 5 (* ' “ T ) ’ 


dx ; 


(4.1) 


where 


^-§ 4 a+^-( c o&a] = 5 4a 5(x-j)6(r-x), 

(4.2) 

Dt gi ° + ax J gj '’ °' 

(4.3) 

^ = 2+c/A. 

(4.4) 
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is the usual convective derivative. Since the solution to this system will depend 
on x x ,y x ,t,r only in the combinations x x - y x andr-r , we put 


(4.5) 

= — r JJ e -^- y ^-^g^ 1 1 y, x)d(t-x)d (x, -y x ), 

(2n) 

where x T = jx 2 x 3 }and y T = {v 2 v 3 } denote the transverse components of xand 
y respectively. The resulting equation forg 4a (x r |j r ;A:,co) has a regular singular 
point at the so called critical layer where 

co-AT/(x r ) = 0, (4.6) 


and the solution has an infinite singularity there, which (as will be shown below) 
produces a much stronger singularity in the propagator y . t (defined below). The 

original singularity in g va (x r |j r ;£,co) can be eliminated by constructing an ‘inner 

solution’ in the vicinity of the critical layer which has continuous pressure and 
normal velocity components but discontinuous tangential velocity components 
(i.e., it has a slip line). However, we do not present the details here because our 
interest is in the corresponding singularity that occurs in the source variable 
j when the observation point xis in the far field. This singularity can only be 
found by working with the 4 th component of the adjoint vector Green’s function 


gav(j> X |*T) = g V a(*T|j,T) (4.7) 

which (unlike the direct Green’s function that is determined by equations in the 
observation variables x,t) is determined by the following set of equations in the 
source variables y,z (Morse and Feshbach (1953), p.870; Goldstein (2006)) 


M 

Dz 


gj4 


dvj 


— c 


8gZ, r-\8e„ 


dy, p 8y: 


S 44 


dg 


54 




= 0 


(4.8) 


M 4 

Dt 


v-Hr-i )g;^ 


dy t 




S(x- y)S(t-r) 


(4.9) 
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(4.10) 


M 4 

Dt 


1^4 

p tyj 


g, 


1 4 


= 0 


where 


d a „ , , a 

= 1- V,. (j) 

Dr ar ab- 


using Appendix A and retaining only terms O(s) leads to 


D o8i4 , _a du ydg44 d g 54 


Qo/4 

Dt 


+ g i 


-c 


aj,. av, ay ; 


= ^ 


U — + U 


V± + V± 


BY By x 1 Svj j 


g 


i 4 


8V i 

-g% — ~gl4 

fyi 


BU (l) „ 8U 


fy t 


+ 5-1 


dY 


+ S n 


J dgl / ~ 1 o a , ^§54 

— §44 + 


a A 


BY R 


BY 


^kJjk = S(x . y)S(t . T ) + e 

Dt By, K K } 


U — + D (1) — + Vj — 
BY 8y l By 


§44 + 


dgt 


a , ^Sl4 


1 J 


BY 


-(r-i)g a 44 


r 8Vj du A 
V tyj + BY j 


Dpgt 4 

Dt 


= £ 


U — + U 
BY 




By j 1 By 


J J 


§54 r g 14 


where j- 2,3 and D 0 /Dt is defined by (4.4). 

4.1 The ‘Outer Solution’ 

Equations (4.12) to (4.14) suggest that g^ should expand like 


(4.11) 


(4.12) 


(4.13) 


(4.14) 
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Cl _ „a,0 +F p.0,1 , 

o v// <5 v// °6v// 1 * * * 


(4.15) 

in any region where it is nonsingular and that the lowest order term can be 
expressed in terms of the single scalar Green’s function G“ (j,r|x,r) by 


a,0 2 

§i4 =~ C 


which satisfies the 3 rd order equation 


d D 0 d£/ d 

\8y t Dx dy t dy l 


a, o 

§44 


DqG° 
Dx 2 


(4.16) 


d ~ 2 d D 0 G a 
v dy< C dy j Dx 


dIg „ a 
— + 2 - 
Dx 3 




r du~ 2 d A 

c G 

dvj dy x j 


= 8(x-.y)8(r-x) 


(4.17) 


in this region. 

Since the coefficients in (4.12) to (4.14) are independent of time and only 
depend on y, through the slow streamwise variable Y its solution should depend 

on x x ,y x ,t and t only in the combinations x, -y x and t - r . And it, therefore, 
makes sense to put 

g" 4 (j r ,7|x r ;A,co,s) 

(4.18) 

= (y^\x,t\s)d ( t-x)d (x x - y, ) 

( 2n ) 

It then follows from (4.16) that 


''a, 0 2- 

Si 4 = ~ C Z 


-((o-Uk)-2k— 

dv i V ' dy i 


G a = 


c 2 dG, 


(to -Ukf dy, 


0 ,i=2,3 


-ike 


g,T =-c-k(a-Uk)G " = 7 —G„ 


( co-Uk ) 




(4.19) 


where 

G X) (y T j\x T \k,ed) = 


12 



(4.20) 


, (co U{y T )k) rr 
( 2 tt) 

= -/(co-t/(j r )/c) 3 G° (j r ,7|x r ;A:,co) 


Is determined by 


4 <4 = 


&(x T ) 

(2k) 2 


4 = 


+ 1 -- 


k 2 c 2 


fyj ( kU - co) tyj ( kU - co) 


y = 2,3 


(4.21) 


(Afsar et al (2006) show that this formulation is computationally advantageous 
when the mean flow is only known at a discrete number of points)Then, 
corresponding to (4.15) 

Svct (jV ,Y\x T ',k,co,£^ = (y T ,Y\x T \k,co) + £g°£(y T ,Y\x T -,k,(D,s) + ... (4.22) 


in the region where is nonsingular, and the leading order term satisfies 


-i(co-Uk ) + ik c 2 g 4 a 4 ° +g a 5 f = 0 


/ ttk\ ««,0 , ~ a .O dU 3 6 sl 4 
ly(X) Uk ) g j4 + g l4 c 

fy t d Yi 


d& afi 

-^- = 0, for/ = 2,3 

fy t 


(4.23) 


-i(a-Vk)gl» 

fyi 


+ ikg a u° = 


( 2 x) 


■ S ( Xj y T ) 


(4.24) 


i(co-Uk)gl f =0 


(4.25) 


there. The fifth component,^; 0 , which is determined independently of the 

remaining components by (4.25), is indeterminate at this stage of the analysis. 

These equations also have a regular singular point in the vicinity of the 
critical layer which passes through the line of singular points, say y T = y c , where 
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(Q-kU(y c ) = 0. 


(4.26) 


The resulting solution will exhibit a strong infinite singularity when the source 
variable ^approaches this point, which can be attributed to the breakdown of 

the parallel flow model. It can be eliminated by constructing a local non-parallel 
flow solution in the vicinity of this point-referred to here as the “inner solution”. 
The critical layer, which disappears when the source variable amoves onto the 

jet centerline, will lie on the level curve/surface of U , say S c , where the dominant 

balance of the outer (parallel flow) solution turns out to be (Wundrow, & 
Goldstein 1994) 


(sTx r ) , . . 

g'a 4 = t + b (T (s,Y\x T ) + c (7 {s,Y\x T )y ± \ny 1 +..., cr = l,t 


a„ , 

i,0 er 


y 


, a0 ^_L4 Y \ X T ) * / v i 

g ±4 = + c 14 {s,Y x r )ln y ± +..., 

y± 


< a ,o a 4 


g 44 


(5,r|x r ) 


y± 


+ b 4 (s,Y\x T )y ± +c 4 (s,Y\x T )y 2 L \ny ± +..., 


(4.27) 


-a, 0 

£54 =■ 


where the a v ,b v ,c v ... are completely determined by the outer solution while the 

dots in the last equation are meant to reflect the indeterminacy of g 5 a ;°at this point 

in the analysis, 5 denotes the transverse distance along S c and y± 

denotes the perpendicular distance to this surface, i.e. , the distance in the 
direction of its unit normal n (see figure 1 ). 

4.2 The ‘Inner Solution’ 

The solution in the ‘inner region’y ± =y ± / 8 = {y T -y c )-n/S = 0(1) , 

where 5 denotes an appropriate scale factor for the critical layer, should therefore 
expand like 


ga 4 =8- 2 g°a4(y±’ S ’ Y ) + d~ l gl4(yL’ S ’ Y ) + g 2 a4(yL’ S ’ Y ) + ---’ C>=\,t 

g a L 4 =^ l gl A (y L ,s,Y) + \ndg\ A (y L ,s,Y) + g 2 LA {v L ,sj) + 5gl A (y L ,s,Y) + ..., 
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gy 4 = 8 ~ l gl 4 (y±’ S ’ Y ) + gy 4 (y±’ S ’ Y ) + Sgy 4 [y L ,S,Y) + . . . , V=4,5 


(4.28) 


Then, since the distinguished limit corresponds to S = 4s , 



1 + 5/ 4 0) y± + \z 2 Kc (s)yl+ — 


(4.29) 


and, for any function / = /( y ± ,s,Y) ,the chain rule implies that 


r 

v 


BY y 


y T 


r W} 

dr 


df dy± _ f df A 


+ 


y±,s 


BY 


8y ± BY 


y±s 


\_dj_ 
5 By ± 


dy c 

BY 



(4.30) 


where y c = y c ft, the subscript c indicates that the quantity is to be evaluated at 
y c and the primes denote differentiation with respect to y ± , the lowest order 
terms in these expansions should satisfy 


Qy± 8y ± 


ikU'y ± g ]4 + ik c; g 44 + g 54 = 


v„ -u 


Sy c 

BY 


d A, 


By 


ikU' c y L g° tA - ±(4 g 44 + g 5 ° 4 ) = 


V„ -U 


BY 




-g<4> 


(4.31) 

(4.32) 

(4.33) 


0 1 


%14 -T- 


1 dh±g t 4 + Bh s g 14 


\K 


Bs 


dy ± 


= 0 


(4.34) 


ikU' c y ± g 0 5 4 = 


V„ -U. 


& 

BY 


a— ^ 54 

Qy± 


u. 


BU 

BY 


+ VU' 


0 


(4.35) 


dh sg[ 4 _ Q 

dy±_ 


(4.36) 
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( 4 . 37 ) 


%14 -T- 


1 Sh ± g] 4 + dh s gl 4 


v h ± ds 


Qy± 


= 0 


ikU' c y L g 0 44 +ikgf 4 —j- 


3 A 


1 dh ± g t4+ dh s g 14 


ds 


dy± 


\ 

f dy ^ 

v -t/ c — 

J 

t c dYJ 


d g i+ dg 


dy. 


dY 


( 4 . 38 ) 


where F c = V(y c ,Y)-n, 




' dy T 

dU 


Js 


• = K M + 5/? l c + • 


(4.39) 


These results were obtained by using (A. 12), (A. 13) and (A. 16) to get (4.35), 
and retaining higher order terms in equations (4.34) and (4.36) to (4.38)in order 
to circumvent the complexity of using the expansions (4.29) and (4.39). 


The first of these equations, which can be written as 


giu: = 


dy 


where 

— ~2 0 , 0 
r = c c g 44 +g 54 , 

can be used to eliminate gj° 4 in the 2 nd equation to obtain 


(4.40) 


(4.41) 


where 


d 2 _ 




y -2a~ 


_ dy _ 


= 0 


(4.42) 


a~ = ik - 


u: 


c 


V. -u 


dy c 

dY 


(4.43) 


The solution to these equations that matches onto the 1 st and last of equations 
(4.27) is 
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( 4 . 44 ) 


and 


00 

y = 2a 2 a 4 c 2 exp.(a 2 v 2 ) J" exp .(-a 2 7/ 2 Jd7j 

y± 


^14 = 


2« 2 c c 2 

U[ 


8 

8 y± 


exp.(a 2 v 2 


00 

) | exp.^-a 2 }] 2 yj] 

y± 


( 4 . 45 ) 


These results, along with those of Appendix B which contains the details of the 
critical layer expansion, show that the inner expansion for the Fourier Transform 

g a vA (y T ,Y\x T ;k,ti),£) of adjoint vector Green’s function g" v can be 

summarized as follows: 


2 * a * a ^CCC / — \ l 1 \ 

A § 44 +g 54 = r w { a y± j + OUs) 

ds v ’ 

( 4 . 46 ) 

gi = 2 

( 4 . 47 ) 

* = 2 Ti" ( “* ) + 0 (l) 

(4.48) 


(4.49) 


§ ±4 


2 a 


yfsikU' 


f \ 8 8c 1 a A 2 ~- ^ 

h, — c -^- + k c:a 


h ds 


8s 


-C 4 


W 


(«Ti) 


(4.50) 


J 


-ik + — 
^ 8Y 




, 8 1 

^14+7" 

h s 


f 


1 dh ± g? 4 + dh s g a ±4 


\ h ± 


ds 


8 y 


( 


4 - J 


v-u. 


8y c 

8Y 


2 a~y ± - — 
d y± V 


r + u' c gl 


= a 


l + --^w(ay ± ) 
a dy L 


(4.51) 
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where equations (4.34) to (4.38) and (B.3) were used to obtain (4.51 ) and we 
have put (Abramowitz & Stegun, 1965, p.297) 

00 

w { a y±) = aexp.(a 2 yl ) J" exp = 

y± 



1 1 

— — + — as ay, -» oo 
2 ay L 4 a y ± 


(4.52) 


and 


a = ikU' c a 4 . (4.53) 

Corresponding to(4.27), G 0 should behave like 


G 0 ~ a(x,Y,s} + b(x,Y,s}y\ +c[x,Y,s^)y\ \ay L +... as y ± — » 0 (4.54) 


where d = G n 


* v ± 


=o 


is defined by (4.53) and 


b=- 


1 d , del a 
-h 


- + k ci 


}l c 2 ds ds 

V n ± U c 


(4.55) 


We are ultimately interested in using these results to obtain a uniformly valid 
expression for the space-time Fourier Transform of propagator (3.5) which, in 
view of(4.7) an (4.18), will depend on g" 4 [y T ,Y\x T \k, (£>,&) and its derivatives with 

respect to y T (see(5.14)-(5.17) below). This will be accomplished by using the 

product rule (i.e., multiplying the inner and outer solutions and dividing the result 
by their common part in the overlap domain (Van Dyke, 1975) to combine the 
inner and outer expansions for g" 4 ( y T ,Y\x T -,k,v ),& ) .The most singular 


contributions come from the derivatives dg“ A / dy } , j = 2,3 ; i * j with a typical 
term being dg a X4 / dy ± . Equations (4.19) and (4.48) show that the inner and outer 
expansions for this quantity are given by -ikd c 2 G 0 / (co - Uk)\ / dy ± and 


g“ 4 =[2ac 2 c a 4 / eSU' c \d 2 w(ay 1 )/ dyl and equations (4.52) and (4.54) show that 

they will match in some overlap domain where y ± <sc 1 and y ± » 1 . It therefore 
follows from the product rule that 
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(0«,V8Vx) unifonnlyvalid =-* 


0 c 2 G 0 
dy ± ( co-Uk ) 



w(ay L ) 


ayl 


(4.56) 


is a uniformly valid approximation to dg a u / dy ± that remains bounded at y ± = 0 

and reduces to the appropriate inner and outer solutions within their respective 
domains of validity. 

5. The Far Field Spectrum 

The focus of this paper is on predicting the acoustic spectrum 


1 oo 

I m (x) = ^-je‘ m ‘p 2 (x, , )d , l (5.1) 

—00 

which is just the Fourier transform of the far-field pressure autocovariance 
defined by 


| 1 

P 2 {x,t o) = — J p\x,t)p'(x,t + t 0 )dt, (5.2) 

where T denotes some large, but finite, time interval. (In most cases, capital 
letters will be used to denote Fourier transforms of the corresponding lower case 
quantities) Then, since p' e ^ p' as x -» oo, equation (3.4) can be inserted 

into(5.2), and the integration variables changed to tj = t - x 2 and x x = t 2 - % , to 
obtain 

P 2 (x,t 0 ) = 


^ i oo 

~yr 1 11 II y vi ( x \y> t + t o~ X K/ { x \yi’ t ~ x i) e vj t) e", (y 1 ,T 2 )dy dy x dx dx 2 dt 

-T V 


+t 0 +x l )y Vil [x\y,,t l )^ vjill {y,y 1 -y,x,)dydy t dt, dx, , (5.3) 

-oo y 


where we have taken the liberty of using the same symbol to denote slightly 
different functions and put 
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y vj (x\y, t -x)^y vj {x,t\y,T), 


(5.4) 


J 1 

Xy/vi [y^) = — \ <■ (y, *o Xi(y + T i» T o + T ) dt o ( 5 - 5 ) 

is the modified fixed frame density-weighted, fourth-order, two-point, time- 
delayed fluctuating velocity correlation and 

n = y,-y (5.6) 

denotes the separation vector between the points. This result can be put into a 
more transparent form by introducing a new propagator y . ,, defined in terms of 

the original propagator y by 

oo 

lyjvi ( x \ y>^ t + x )= J +^ + t )y^(*|t + Vi)<*i (5.7) 


to obtain the following, more compact, expression 

oo 

A*.')- 1 II Yv jv i (*| y;tj,t+r)^ vM (y;tj,r)dy dy dr (5.8) 

—oo V 

for the far-field pressure autocovariance in terms of the four-dimensional two- 
point, time-delayed momentum/enthalpy flux correlation tensor K,Ar.n,t) ■ 

This result is essentially exact, but the source correlation can contain purely 
convective, and therefore non-radiating, components. In an effort to minimize 
such components (in the Lighthill equation context), Ffowcs Williams (1963) 
introduced a moving frame correlation tensor which, with the present formulation, 
is appropriately defined by 


Km ( t ; r) = z vml ( y; % + \u c z, t), 


(5.9) 


where 

S = U C T (5.10) 

denotes a moving frame coordinate system, with U c being the “convection 

velocity” of the turbulence. And since this objective is consistent with the 
philosophy of this paper (as outlined in the Introduction), we insert (5.9)into (5.8) 
to obtain 
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(5.11) 


P 2 ( x ,t)= J |{y v j l u(x\y;$ + iU c T, t +T)%“ vJ (y;$,T)dy d£ dr . 

—oo V 

The modified fixed-frame correlation (j;t],r)is related to the autocovariance 
of the generalized velocity/enthalpy Reynolds stress tensor 


I t 

Kjvj (t; t) = — J [pv>; - pv>; ] ( y, z 0 ) 


P'V’/ - P'V ; / ] (y + %To + T ) dz o (5.12) 


(which is about as close as you can get to what is actually measured in turbulent 
flows) by the simple linear transform 


*vj\u ~ %) xl + 5 yil R vjkk ^ + ^ Y 


^y/^| xl^iikk 


(5.13) 


with an identical relation holding between the corresponding moving frame 
correlation tensors ^ j[yl {y,^,z) and R^^(y;^,z). 

The far-field spectrum can be calculated by taking the Fourier transform in 
time of (5.1 1 ) and noting that (5.7) is of convolution form (Morse and Feshbach, 
1953, p.465) to obtain 


7 «»(*W = 27t J { r y / (^|T;a>) r * / (x|j + ^ + /C/ c T;(o)e mT ^,(y,^,z)d^dz (5.14) 

—oo V 

where 


r,- J e‘^ , ~ z \ v j(x\y,l-z)d(t-x) (5.15) 

—00 

is the Fourier transform of y . and we have introduced I m (* | j) ,the acoustic 
spectrum at xdue to a unit volume of turbulence at y , i.e. , 


/ co(-* : ) = | / a ) (^|T)^ ( 5 - 16 ) 

v 

in order to make the formula more transparent. 

Appendix A and equations(3.5),(4.7),(4.18) and(5.15) imply that 


r 

vy 


oo 


( x \ y’<°)= | 

—00 ■ 


_°_ e ik i x i-xi) 

dy t 


gU{yT’ Y \ x T; k ’®’Z) 
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dk 


(5.17) 


- (y - 1) S vl /' (xi - Vl 1 g« 4 (y T ,Y\x T ; k, to, e) 

fyj J 

to lowest order of approximation, where the integration contour would have to be 
deformed below the real axis in order to pick up the contribution from the 
instability waves if the Green’s function were required to be causal (Goldstein & 
Leib, 2005). But our most recent studies suggest that this contribution should be 
fairly small for cold (i.e., unheated) jets at the relatively low supersonic Mach 
numbers being considered here. Then, since our focus is on jets of this type, we 
tentatively neglect it by requiring that the integration contour lie along the real 
axis. It can easily be incorporated, especially if it is assumed that it is 
uncorrelated with the non-causal contribution. 

The result is much simpler in the far field where U (x r ) -» 0, c 2 -» = 

constant, and 

e -x r ^k 2 -((o/c 00 ') 2 

§v 4 ^ / == #v4 (cp,^,co| y T j) as x T — » oo, (5.18) 

Jx T 

where we have used (4.5),(4.7),(4.18)and (4.21) (with y , replaced by 

Xj), x T = |x r | and cp = tan 1 x 2 / x 3 denotes the circumferential angle. Inserting this 

into equation(5.17), and using stationary phase to evaluate the integrals, shows 
that 

r v,(*kl. y ) = 



2ra'cosin9 d 


#v 4 — cos 0, co y T ,Y 


-(y-l)Sy, dU [ r,fT> ^ t ( <p, — cos 0, ra |/ r , rjl e - '"" 1 

V J 


where x = |x| , 0 = sin 1 (x T lx) denotes the polar angle measured from the 

downstream jet axis (see figure 2) and the term in curly brackets should be 
interpreted as an operator that operates on everything to its right. So that, in 

particular, d/dy f acts one>“ iO, - VlCOS0/Cx ^ 4 . This result is valid for all values of y T , but 
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when y T does not lie within the critical layer, equations(4.19) and (4.20) can be 
inserted to show that 


r ij( x \yT> Y ) = \e 


iv\ cos 0co/ c rfl d c 2 /' co 2 d cos 9/coo 


iy-y cos Ucd/ Cqq 


tyj [M ( , 7) cos 0 - 1] 2 


i8 v (y-l)dU 1 

co dry j [M(j r ,7)cos9-l] 




co 


cp, — cos 9, co y T ,Y 


V c « 


(5.20) 


F 4 • (* I j r , f) = —e ib>yi Cos0 / Cco — 

7 v ’ co dvj 


where 


3 -fcoy, cos9/c x 


^0 


CO 


cp, — cos 9,(0 y T ,Y 


V L O0 


M(jj,,y)cos9 -1 


M(y T ,Y) = U(y T ,Y)/c x 

denotes the “acoustic” Mach number, (cp,A,co|j r ,F)is defined by 


(5.21) 


G 0 -» 


-x ± ^k 2 -((o/c x ) 2 


^ 0 (cp,/c,co|j r ,7), asx ± H>co 


(5.22) 


and the term in curly brackets should again be interpreted as an operator that 
operates on everything to its right. 

This result is important because it provides an explicit representation of the 
singularity that will occur in the propagator when the parallel flow model is 
invoked. It shows, for example, that there is a singularity in the transverse source 
variable ^whenever 


M(y T ,Y) = 1/ cos 9. 


(5.23) 


and that the strongest singularity inr ;> , which results from differentiating 

the first term in(5.20), isO(v^ 3 ). This is a strong non-integrable singularity that 

will cause the far field pressure to become infinite and must, therefore, be 
eliminated (by using the uniformly valid non-parallel flow components given in 
Appendix C) before meaningful predictions can be made. 

Streamwise correlation lengths seem to be very long in certain regions of the 
flow (J. E. Bridges, personal communication) and it is therefore necessary to 


23 



account for retarded time variations in that direction. But the transverse variation 
can also be important at sufficiently high Mach numbers. We assume that the 
transverse correlation lengths are small compared to the cross stream 
dimension, say D, of the mean flow. ThenT*, (x|j r +g T ,Y; to) , which scales with 

D and the acoustic wave length c x / co , can only vary significantly over the 

correlation volume when co Dt is large. In which case F* ; (x|j r +g T ,Y; co) can 

be represented by its high frequency or WKBJ approximation* (Khavaran, 
Bridges & Georgiadis, 2005) 




/— S(<p,0|^,r) 


(5.24) 

where (x|j r ,7;co) expands like 


^J ) (x|j r ,7) + — ^ 1 / ) (x|j r ,7)+ — A [ ^(x\y T j) 


“44 




4 2 ), 


V c <» y 


(5.25) 


and S(cp,0|>v,7), which satisfies the Eikonal equation 


(Wy T s)-(v yT s) = [l - M ( ) cos e]" c 2 Jc 2 - cos 2 0, (5.26) 

varies on the scale of D. It is, therefore, relatively constant over the correlation 
volume and can be expanded in a Taylor series for variations on this scale to 
obtain 


r iv (*|.Vr +£ r ,T;co) « r* ; (x|j r ,F;co)exp 


CO 


^•V^5(cp,0| Jr ,7) 


(5.27) 


which allows us to write the formula (5.14) for the far field spectrum in the purely 
algebraic form 

7 »(*W 


-» 


*Actually, the solution can be a sum of terms of this form, but the final result turns 
out to be the same and we therefore consider only the single term for the sake of 
simplicity. 
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' 2n Y 
V X ) 


2na> 


siner vy (x|jv)r;^x|j; r )0* vAl/ j;-^-cos0, — V y S,co(l-M c cos0) 


(0 

C„ 


asx — >oo (5.28) 

which depends on the turbulent source correlations only through the spectral 
tensor 


Kj»i{y; k 



i(kfy+k T -^T) —M 
e />/ 


(y,£,i)d£dx. 


(5.29) 


where the asterisk is being used to denote complex conjugates. The comment 
following equation (5.23) implies that the strongest critical layer singularity in this 
result is now O^y ] 6 ) when non-parallel flow effects are neglected\ 

The presence of the enthalpy source e ' Aj , significantly increases the complexity 

of these results. But equations (2. 13), (2. 14) and (2.23) suggest that the 
dimensionless ratio, e' 4j /U J e' j (which determines the relative importance of 

e' 4j and e' tJ in the acoustic analogy equations, where Uj denotes a characteristic 

jet velocity)should be 0{v[IU ,) for cold, i.e. unheated, jets when the Mach 

number is 0 (\) , (see, Lilley, 1 996 and Morfey, Szewczyk, & Tester, 1978), 


because c 2 ' is expected to be 0(v/ ) 2 in this case. The enthalpy source e' Aj should 

then be small and can, therefore, be set to zero if we focus on cold jets — which 
we now do for the reminder of the paper. It is, of course, possible that equation 
(5.28) can still produce reasonable predictions even for hot jets. But this would 
require modeling many cross-coupling terms that would be difficult to estimate 
from any standard RANS based code. It may, however, be possible to do this by 
extending a RANS code to yield information about the RMS temperature 
fluctuations (as was done by Khavaran and Kenzakowski, 2007) if the velocity 
and sound speed fluctuations are assumed to be uncorrelated in the sense that 


0> : 




y;— cosG, —V yT S,Gi(\-M c cos0) 

Crr Cm 


= 0 


(5.30) 


when p * v and either p or v=4 -but the result may be worse than neglecting 
e' 4j entirely, because the velocity and enthalpy (sound speed) fluctuations may 
interfere destructively (Freund, 2002).We, therefore, assume that e' 4j = 0, i.e., that 
(5.30) holds whenever p orv is equal to 4. Then since (5.12) and (5.13) show 
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that 0* ijkl is symmetric in its first two and last two indices, equation (5.28) can 
now be written as 

4(*| y)^> 


'2n\ 2 
V x ) 


sin 6 G tj (x\y T )G* kl (x|j r )<D * ljkl 


y;— cos0, —VS , co(l -M c cos0) 

v Ao C co j 


where the symmetric propagator 


asx — »oo (5.31) 


°ij ( x \yr ) = Gji ( x \Y t ) = 2 [ T v (* \ y T ) + r Jt (* \y* )_ 


(5.32) 


still has six independent components, but only certain combinations of these 
appear in the simplified formulas (see section 6.1 , below) used in the actual jet 
nose predictions presented in section 8 of this paper. Appendix C uses the so- 
called “product rule” (Van Dyke, 1975) along with the results of Section 4 to obtain 
appropriate uniformly valid expressions for these quantities. 


The critical layer lies at the transverse location y T = y c determined by(5.23). 

Figure 3 is a sketch of its overall configuration for a round jet. It shows that the 
critical layer usually moves inboard from the lip line toward the jet axis with 
increasing downstream distance and fixed 0 and that, while the initial migration 
tends to be relatively slow, it quickly moves onto the jet centerline and 
disappears just after reaching the end of the potential core. The right angle 
intersection is shown because U = U(y T ,Y) and 


U' c = dU / dy 1 \ y ^ =yc =0 for y c = 0 (5.33) 

in this case while the data seems to suggest that (dU / d 7) = 0(l) . So it follows 

from equation (B.2) that dyj dYrc\us\ be infinite there. This behavior is expected 

to be generic because the on-axis mean-flow will usually behave as if it were 
axysymmetric. 

The upstream critical layer thickness, say5 f , is also shown in figure 3. 

Equations (4.43) and (4.52) imply that it will become infinite when the critical 
layer reaches the jet axis, which would cause the critical layer solution (4.46) - 
(4.51 ) to completely cancel out the outer solution singularity there. But equation 
(4.32) shows that the large value of dyj dY causes the critical layer expansion 

(4.28) to break down before this can occur and a new asymptotic solution will 
apply in the so called blending region indicated by the circle in figure 3. But 
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actually constructing this solution would be very difficult because the relevant 
flow structure (including that of the mean flow) is no longer quasi-parallel. 

The weakly non-parallel flow structure appears to reemerge further 
downstream. But the outer parallel flow solution will remain invalid in the vicinity 
of the jet axis as long as the difference, sayA(7 , between the centerline velocity 
and the critical layer velocity U ^ ^ = c x / cos# remains of order of the jet spread 

rates. We refer to this inner region as the modulated layer. The asymptotic 
structure and scaling of the relevant Green’s function solution will differ from the 
critical layer structure in this region, but we do not derive the result here. The 
detailed analysis is similar to, but much more tedious than, that of section 4. The 
primary difference is that the inner solution is now determined by a second order 
partial differential equation rather than by the ordinary differential equation(4.42). 
Fortunately, this equation can still be solved in closed form, but the solution is 
considerably more complex than the previous result. We can, however, simplify it 
by noticing that the streamwise length of this downstream region is typically 
much shorter than the streamwise length scale of the mean flow and that the 
RANS solutions show that the mean flow stagnation enthalpy H 0 =c f + U 2 / 2 is 
very nearly constant. The end result is that y turns out to be proportional to 


Yo = <^ X °^ q r(x 0 -l,q) (5.34) 

where r(x 0 — 1, q) denotes the incomplete Gamma function with argument q 
given by 


ikJJ” v. , , cocosO 

_ 0 and r _ 

4F ° c 


and scaled modulation parameter 


Jo 


-* 0 (aC/ + £/('>) 
2F 


where 


AU = \U.-U , \le. 

cr. lay. 


(5.35) 


(5.36) 


(5.37) 


The appropriate lowest-order expressions for the relevant propagator 
components can be obtained by making the substitutions 

A) ^ - X 0 (5-38) 
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YofeXo) 


(5.39) 


w 


'( z o) 


1 f a 1 (y-i)5 du c 1 

2|dq 2K,'cn 57 q 


(i-Xo) 


_a _a_ 

5x 0 5q 


in equations (6.27) (below) and (C.18) to (C.22), where z 0 is defined by (6.15) 
(below). This solution merges smoothly into the parallel flow solution in the 
downstream potion of this region where |x 0 | becomes large. But it must also 

merge onto the blending layer solution at the upstream part of the region. We do 
not attempt to construct this solution here -which cannot, in any case be done 
analytically. But since the thickness of the upstream critical layer is infinite when 
it becomes perpendicular to the streamwise axis, we expect the blending-layer 
solution to also spread out into the axial direction and smoothly merge into the 
downstream region also shown in figure 3. 

6. Source Modeling 


The previous sections were primarily concerned with the purely acoustic 
effects, frequently — but somewhat misleadingly- referred to as propagation 
effects. The remaining ingredient involves the (near field) fluid mechanical 
effects, which are inputted by modeling the turbulent sources that enter the 
acoustic formulas through the modified spectral 


tensorO’ 


ijkl 


f \ 

y ; — cosG, — V >v 5',co(l-M c cos0) 

^ ^ " . 


Ideally, we would like to model 0* A/ directly, but the models must be based 


on experimental data and the experimentalists are unlikely to provide the 
requisite data in the near future. A possible alternative is to develop models for 
the ordinary spectral tensor 


' ¥ ffld(y; k i> k x* 


co) = — f 

’ 2n J 


V 


\kfc\+k T '%T) 


Rf kl {y,Z,x)d£dx, 


( 6 . 1 ) 


which is implicitly related to 0* ijkl by (5.13). But even this quantity has almost 

never been measured (see, for example, Harper-Bourne (2003))-and then only 
at very low Mach numbers. Various components of the two point correlation 
tensor appear to be the most extensively measured turbulence quantities in high 
speed air jets. We are, therefore, forced to develop models for these quantities 
and calculate the spectra from(6.1 ). 

There are two main requirements for such models. The first is that the large 
number of independent components (there are 45 of these in all) be brought 
down to a manageable level and the second is that realistic, but still relatively 
simple models, can be obtained for those components. The first of these (which 
is discussed in the next section) is usually accomplished by introducing add hoc 
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symmetry and/or statistical assumptions such as local isotropy and quasi- 
normality. But, as indicated in the introduction, such assumptions are not viable 
at the high Mach numbers of technological interest. A more realistic approach is 
therefore introduced in the next section. Specific models for the remaining 
spectral components are introduced in section 6.2. 

6.1 Kinematic Modeling 

As indicated in the introduction, the usual quasi-normal and local isotropy 
assumptions are untenable at high Mach numbers, but developing a viable 
formulation for a general non-normal axisymmetric source model 
(Batchelor, 1953; Goldstein, and Rosenbaum, 1973 b, Kerschen,1983) is 
enormously complex and the result would probably be too complicated to be of 
practical utility. We therefore proceed by developing a non-normal but isotropic 
model and an axisymmetric but quasinormal model. 

6.1a General (non quasinormal) Isotropic Turbulence 

Equation(5.13) shows that ^ will be an isotropic tensor if Z?" is an 

isotropic tensor. Batchelor (1953) points out that the most general fourth order 
isotropic tensor is of the form 


+ BUA, + CU* 6, + DU, 5 , + E^L k 5, + F^L, 5, 


( 6 . 2 ) 


+ GUA J + Hd iJ d kl +Id ik 8 ] i+Jdi l djk 

where A,B,...J are functions of y , t and |£| . But it follows from the definition of 
that 


and, therefore, that 





(6.3) 


K, = jZAAA + buAi + + G ^A% + (5,s, + 5,s Jf ) 

+ (e. t 6„ + ak ) + n\b u + /(8„s„ + s„s Jt ) 


(6.4) 


which can be integrated over the unit sphere to show that 


J (y;S ,f) dt = J,8„8„ + 1, (S„8„ + 8„8 J( ) (6.5) 

V 

where 
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( 6 . 6 ) 


I x (y;t")=An\ 


0 L 


-^A^+^(B + G)^+H 


edt 


=J<, df-2l*fi a (y;f,r) a 

V V 

and 


( y;?) = 4tcJ 


^At;‘+^(C+E)t; 2 +I i; 2 <«; = 




Inserting this into equation (5.29) shows that 

0* w ( = <i>i 5 A/ + k (8,S /7 + 5 ;/ 5 y7fc ) 

where 


$1 = 0* im (j;O,0,co)-2®* 1212 (j;O,0,co) 


$2 =°*2i2 (j;O,0,co) 

with 

O *1212 ( j;O,0,co) = O * 1313 ( j;O,0,co) 

It is reasonable to assume that d>*.^ (y;k l ,k T ,(o) has the same symmetry, 
that 


where 


and 


0* H (y;k v k T , co) = ^S..S, ; + <j) 2 (5, ,8, + 8„8 jk ) 

^1 =®*111 (y;K’k T ,(0)-2O* m2 (y;k l ,k T ,(o) 
^2 = 0 *1212 [ y ; K ’ k T , co) 

0*212 (y;K>k T > (°) = o* 3 i 3 {y;k u k T f co) 


(6.7) 

( 6 . 8 ) 

(6.9) 

( 6 . 10 ) 
.e. 

( 6 . 11 ) 

( 6 . 12 ) 

(6.13) 
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Inserting these results into equation(5.31 ) and using Appendix D (which contains 
the relevant details of the analysis) as well as (C.18) to (C.22) yields 
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(6.14) 


where the primes denote differentiation with respect to 

z 0 =a 0 y ± , (6.15) 


a 2 denotes the parameter a 2 (defined by (4.43)) with k = cocos 0/c^, C c is given by 


c„ = 


-C 2 2zq>v'(z 0 ) -c 2 2zq>v'(z 0 )cos 2 9 


l + w’ (z 0 ) c 2 ( Jy ) cos 0 — l] 2 [l + w'(z 0 ) 


(6.16) 


(with C defined by (C.4) and the replacements (5.38) and (5.39) for z 0 and , respectively 
inserted after the critical layer disappears) and 


(T) * = (I) * 

1 mi - 1 mi 


y;— cos0, — V y T S ,co(l-M c cos0) 
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V C 00 


(6.17) 


are two independent components of the Doppler shifted turbulence spectrum that 
can, in principle, have entirely different functional forms. 


6.1b Quasinormal Axisymmetric Model 
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We now put 


(y,4,^) = ^=\ Jpv' i (y,t)-Jpv' J (y + £ + iU c T,t + T ] )dt (6.18) 


and assume that the turbulence is quasi-normal(see Batchelor, 1953; 
Goldstein, 2004), which means that 


B*(y,{,T) = Rg(y;{,T)R*(y;f,T) + R* {y^,x)R%{y^,x) 


(6.19) 


since (Goldstein & Leib,2005) 

1 T 

— J pv>) (y, H W v /' (y + n^o + T )dr 0 


= R m (y» n, t) + Ry (y; o, o)r u (y + % o, o) 


( 6 . 20 ) 


We now require the turbulence to be axisymmetric (Batchelor, 1953) in the sense 
that 


K 0^) = A tej + BS ff + C8 u 8i j +D{ 8^, + S ly ^ ) 


( 6 . 21 ) 


where A, B, C, and/) denote functions ofj, z , ^ + £3 and 4ii 

4, B and C are even functions of the latter quantity while D is an odd function. 
This model is chosen because it is the most general of those whose 
mathematical properties have been studied in the literature and because it 
reflects the fact that the cross flow velocity components tend to be much more 
similar to one another than to the stream-wise component (even for non- 
axisymmetric flows). 

It follows from equations(5.13), (6.19) and (6.21) that 
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( 6 . 22 ) 
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where 




(6.23) 


while appendix E (which contains the relevant details for this model) and 
equations(C.18) to (C.22) imply that 


W /*,7 


® \C,f + 2® *. m »(C c (l - C ' c )) +|1 - C, f ® * K2J ] | <?„ [1 + W (z„ )] 


+ 2 ^ ^ 22^22 ( 7 2 22 J 1 2222 ) W ( Z 0 ) + 


(6.24) 
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where C f is defined by equation(6.16). As in the preceding sub-section, we now 
assume that this result along with the symmetry conditions 


(T) * = 0 * <T) * = 0 * <T) * = 0 * 

^ 2222 ^ 3333 ^ ^1122 ^ 1133 ' ^ 1212 ^ 1313 


(6.25) 


also holds when J ^ kl {y;^,r)d^ is replaced by 

V 

Jexpi-^(^c 0 se + § r -V yT 5)^(^;§,T)^. 

V ^00 

Then, inserting (6. 24), along with equations(6.19),(C.18),(E.3) and (E.4) into 
(6.22), and using the result in equations (5.31), (5.29) and (5.13), shows that the 
far-field acoustic spectrum 7 ro (jc|j)can now be expressed in terms of five distinct 
components of the spectral tensor 


® *ijkl - ® *ijkl 


y;— cos0, —V } , t S, co(l -M c cos0) 
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(6.26) 
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by 


-1 — sinfl [®* lm |C 0 | 2 +2® *„, 2 ®(c, (l-C c *))+® * 2222 |l-C c | 2 te[l + v/(z 0 )]| 2 
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+ ( I ^23 I &r 22 r 33 )z 0 >v(z 0 ) ( < &*2222 ® *2233 ) 


with the replacements (5.38) and (5.39) being made just downstream of the critical 
layer. And since 

^ *1111 = ® *2222 = 2 ® * 1212 +® * 1122 and O * 2233 = O * 1122 ( 6 . 28 ) 

for isotropic turbulence, this result reduces to equation(6.14) in that case, even 
though the latter applies when the turbulence is not quasi-normal. This shows 
that the quasi-normal hypothesis does not effect the isotropic limit of equation 
(6.27). An alternate representation of the source terms was given by Musifar 
(1993,2006) 

Equation (6.27) can be rewritten by adding 4 iRe C c (l-C* j 0* 1212 to the first 

term in square brackets and subtracting the corresponding result from the 
second. The modified first term will then be proportional to 

0* mi |C c | 2 +2»rC c (l-C*)](20* 1212 +0* 1122 ) +|1-C c | 2 ®* 2222 , which reduces 

to ®* 11U when the turbulence is isotropic, and in the more general non-isotropic 

case (see Equation(6.16)), becomes identical to the self noise term found by 
Goldstein & Rosenbaum (1973 b) in the Lighthill equation context when the 
acoustic Mach number M is set equal to zero (see Ribner (1969)). The remaining 
(modified) terms can be shown to vanish in the absence of mean shear (i.e., 


(6.27) 
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when M =0, c 2 = c 2 and ^7 0 is determined from the free space Green’s function). 


The vanishing of the last term is a direct consequence of replacing (7 0 by the free 
space Green’s function while the second term vanishes because the subtracted 


C c (l-C c *)' 




term cancels with 


(^ 0 /co) 


in this case. It would, 


therefore, not be unreasonable to follow Khavaran & Bridges (2004) and refer to 
the remaining (modified) terms as the shear noise terms. Notice that they are 
both proportional to<J>* 1212 when the turbulence is isotropic (but not necessarily 


quasi-normal, see(6.27)). This latter spectral component becomes proportional to 
cp * 1U1 when the turbulence is also assumed to be quasi-normal, but cp * 1212 and 


cp * mi can be independently specified in the more general case being 

considered here (which implies that the shear and self noise spectra can be 
independently specified). A similar decomposition was made by Musifar (1992). 
The 1 st term in equation (6.27) can also be rewritten by collecting coefficients 

of |C c | 2 , lReC c and 1 . It can be shown that it is the same as the result that would 

be obtained if the modified tensor were assumed to be quasi-normal with 

the corresponding 2 nd order tensors being axisymmetric. 

Equation (6.27) contains five modified spectral tensor components. 
Equations(5.13), (5.29) and (6.23) imply that they can be expressed in terms of 
five distinct components of the spectrum (6.1 ) of the Reynolds Stress 
autocovariance tensor, which, with the present methodology, can be specified 
independently of one another. It follows from (6.16)and the far field behavior the 
Fourier transformed Green’s function, that the coefficient of O * 2222 and the 

second term in the absolute value squared multiplying <P* 1212 are the only two 

terms in (6.27) that do not have a cos# or a cos 2 0 factor. The entire 90° sound 
field is therefore produced by these terms. But our numerical computations 
indicate that the contribution from the latter term is, as expected, negligibly small. 


As written equation (6.27) is the sum of three distinct terms: 1 ) the term 
containing the square bracket (which we refer to as term / ), 2) the term 
involving ®* 212 (which we refer to as term II), and 3) the remaining term(referred 

to as term III ). The first of these only involves the 1 1 1 1 , the 2222. and the 1 1 22 
components of the Reynolds stress autocovariance tensor. The second only 
involves the 1212 component and the last the 2222 and the 2233 components. 
Terms / and III , therefore only depend on the covariances of the squares of 
velocity components, whereas term //only involves the autocovariance of the 
cross velocity components. Notice that the first member of the latter term has 6 
powers of the Doppler factor 1-Mcos0 in its denominator when the implied 
differentiation is carried out — which will cause it to produce a very large 
contribution to the small angle sound field at supersonic Mach numbers. (The 
other terms in the equation only have at most 4 powers of the Doppler factor in 
their denominators.) 
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6.2 Spectral Modeling 


Experiments (Bridges & Podboy, 1999; Harper-Bourne, 2003) suggest that R^ kl 

cannot be expressed as a simple product of its space and time components ( as 
was often assumed in the past, e.g., Goldstein and Rosenbaum, 1973 b) but that 
a good local model, which is consistent with axi-symmetry but does not account 
for the potentially significant long range coherence effects that are most effected 
by the streamwise retarded time variations, might be 

«,"«= ^ X+? (6.29) 

where 

= /|, 2 +/J 2 with | 2 S |, 2 +f 2 , 

(6.30) 

/ 3 2 = /? 2 (i r ) and y 1 = y 2 (| 7 . ) are arbitrary functions of | r which vanish at 
| r =0, a 00 is a parameter , the tildes denote the normalized variables 

|j = ^ / l x , / l T , t = tA , with = yji; 2 + and the intermediate length and 

time scales /, , / ± , and 1 /A, can be related to the integral scales measured by 
(Bridges & Podboy, 1999). 

Notice that R^ kl (j,£,,0,O) and R^ kl (y, 0,0, r) exhibit the usual cusps at 

£,j=0 and x = 0, respectively, that are observed in virtually all high Reynolds 

number experiments (Pope, 2000, pp. 67, 70, 71,110, figures 3.20, 3.22 (a), 
5.13). Of course, these correlations must, in reality, be smooth functions of their 
arguments even at these points, but the corresponding radii of curvature turn out 
to be of the order of the Taylor microscale (Pope, 2000, pp. 198-200) which would 
certainly be negligibly small at the high Reynolds numbers that are of interest 
here. 

Since (5.13), (5.29) and the symmetry conditions (6.25) imply that 

R^ 222 (y,^ T ,r)caa be replaced by i [ R^ 222 (y,^^ r ,x) + R ^ 22 (y,^, , t)] in 

O *2222 k T ,co(l-M c c°s 0 )), \hatRff 22 (y,^,^ T ,T) can be replaced by 

\ R \\22{y^V^T’ X ) + R \ I m{y^\’^T>' 1 ) in 0*1122 COS0 ))’ etC ’ 

and since these latter quantities depend on only through % T = ^ + E,] , a 

reasonable global model that has the same local characteristics as (6.29) but is 
less compact (and, therefore, able to reproduce the negative loops in the 
correlation tensor (Bridges & Podboy (1999)) would, then be 
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(6.31) 


where the a mJ are parameters, the operators D n K , k = r,l which are defined by 


A - li D r - r ^ , D° = 1, Z>” = D k ...D k , n times, (6.32) 

d£i or 

satisfy the commutation relations D K D A = D A D K . 

Then since (5.12)implies that 
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(6.33) 


must be different for the different components of R^ kl it is reasonable to expect 
that the remaining a m , will be different as well. 

The corresponding spectral tensor is 


^ ijkl 
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X a m l D’" D T e~ x+r d^ T drdt, x 
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(6.34) 


which, as shown in appendix F (which contains the details of the derivation), can 
be written as 
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(6.35) 


where the operators D, , D a are defined by 
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we have put k x =k x l j, a = ml A, k T =k T l T 




2 - 2 
i + © , 
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(6.37) 


(6.38) 


. This shows that the 


is just the Hankel transform of (l + R 2 ) ' /_ exp ~/3\Jl + R 2 + y 
acoustic spectrum (5.31) depends on V 5 only through the square root of 

cos6] c 2 /c 2 - cos 2 0, which can be negative for 

sufficiently small cos0 . 

When |3 and y are set equal to 


P 1 = (lr f + ff |2 ' s) (l r f so. f = « |2 ' s) (lr f s 0 


(6.39) 


where 0 < 5 < 2 and a are real parameters, the span-wise component of (6.29) 
behaves like 
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as \ T — > 0 (6.40) 


and equation (6.38) becomes 
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where 


H(A,B,Z) = (l + B 2 f ' 2 ^J 0 {Brj)e ^ 2+z ^ V+ifV ^ ( 6 . 42 ) 

0 

Z = a(l + i? 2 )' 2 , (6.43) 

(6.44) 

Gradshteyn & Ryzhik (p.7 1 2, # 6.623, 2 (1965)) show that 

h(^z ( 1s) I (l + R 2 ^,B,Z^ = 1 , when Z = 0 in which case (6.35) simplifies to 


\ 67T/i irp 

> (0 ) = —r- 


, (-1) 


" + ‘ D m D l 

U k t U 0) 


(l + R 2 + k 2 r j 


Unfortunately, equation (6.42) can be very difficult to evaluate in the general case 
where Z ^ 0 , but since a will always be very small, of the order of 0.3, 

Z (2 ~ s) /(l + R 2 ) = a (2 /(l + R 2 will certainly be small as well and 

h(z (2 ~ s) /(l + i? 2 ),5,zj can be approximated by 

H(Z (1 ~ S) /(l + i? 2 ),5,z) «H(0,B,Z) (6.46) 


1 2 /s — 1 ] 

Changing the integration variable to Z l ~ q in (6.42) shows that H(0,B,Z ) 
becomes independent of B when Z becomes large. Then since H{Q,B,Z) is 
independent of its second argument for both small and large values ofZ , and 
since k T , which is proportional to the turbulence Mach number, is usually fairly 

small compared to R , it seems reasonable to approximate h{z [ 2 ~ s) /(l + i? 2 ),5,zj 
by 



(6.47) 


which depends on the single argument Z and is much easier to evaluate. 
Equation (6.35) then becomes 
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When s = 1, equation (6.47) can be put into the form 


H(0,0,Z)=\ 




(6.51) 


which can be expressed in terms of the Struve and Bessel functions 
H x (z) and N x (z) , respectively, since (Gradshteyn & Ryzhik (1965), p.31 6, 
#3.366, 3) 
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(6.52) 


and it follows from (Abramowitz & Stegun, 1965, pp.361, 

496,497,#s9. 1 .30, 1 2. 1 . 1 3, 12.1 .29) that all of the higher order terms in (6.48) can 
be expressed as higher order Struve and Bessel 

functions H n (zj and N n (z) , n = 2,3,... plus simple polynomials, which means 
that it is essentially given in closed form for this case. But the result can be 
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rapidly computed even in the more general case where s is not equal to 1 , by 
simply tabulating H(0,0,Z) against its single argument Z. 

All of our computations were carried out with s set equal to 1 , in which case 
(6.40) implies that the transverse correlation tensor R^^y, 0,^,0) behaves like 

1 /- n 3/2 

exp.-(^ r /a j and, therefore, has zero slope but large curvature at this point. 

Figure 4 shows the predicted density-weighted fourth order correlation 
R uu{y> 0,4 r ,0) based on the functional form (6.39) with the s, a values of 1 ,0.1 

that were used in the actual acoustic computations. It also shows the 
corresponding quantity inferred from the Bridges and Podboy (1999) 
measurements of the second-order correlations by using a locally homogeneous, 
quasi-normal approximation. The measuring point was on the nozzle lip line 7.5 
diameters downstream of the exit in a Mj =0.9 jet. 


7. Application to a Round Jet 


Current aircraft are usually configured with circular engine exhaust streams. 
We, therefore, make the present results more concrete by considering the case 
of an axisymmetric mean flow, where r ± = |y r | , cos cp = y 2 / r ± , 

M =M(r ± ), c 2 =c 2 (r 1 ). The coefficients in equations (4.21) are expressed in 
polar coordinates in Appendix G and the linear operator in (4.21 ) now becomes 
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which implies that g o will be of the form 
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It follows that 


(7.2) 


(7.3) 


(7.4) 
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for r ± > r[ (and a similar equation with r ± and r[ interchanged when r ± < r' ), where 
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of wjy and wf * , is independent of r ± . It follows from (4.17) and (4.21 ) that 
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Since the spectral functions are expected to be independent of cp' in 
equation(6.27), it is appropriate to average it over this angle. Appendix G and 
equations (C.18) to (C.22) can be used to show that the result is given by 
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with C c given by (6.16), 
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(7.14) 
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and 


D = co(l-Af cos0) 


(7.18) 


where the Doppler weighted cosine function C c is defined by (6.16) with C given 

by (C.4) and the replacements (5.38) and (5.39) inserted after the critical layer 
disappears. 

8. Prediction of the Sound Field and Comparison with 
Measurements 


In order to use the point-wise acoustic spectra (7.12)-(7.17) to predict the 
radiated sound field it is necessary to relate the modified spectra 

f - \ 
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to the ordinary moving frame spectra 


by inserting (5.13) into the basic 


definitions (5.29) and (6.1). Appropriate special cases of the source model(6.47) 
and (6.48) can then be inserted into the result. 

The simplest formulas that are consistent with the Bridges & Podboy (1999) 
measurements correspond to setting the integer 5 equal to 1 in (6.47) and 
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retaining only the lowest-order source non-compactness term, i.e., the a 10 term. . 
Then since 


equation (6.48) reduces to 


^ d _ k x a d 
1 8k, ~ Z 8Z 


( 8 . 1 ) 




( 1 \ 3/2 

Z\\ + R 2 +k 2 ) 


I ( a 0,0 «l,o) 


3//(0,0,Z) a 2 dH (0,0, Z) 


l + i? 2 +£ 2 Z 


dZ 


( 8 . 2 ) 


+« 1 , 0 ^ 1 2 


15/7(0, 0,Z) 
(l + i? 2 + k 2 ) 2 


a 


■ + ■ 


6cr 


Z J Z(\ + R 2 +k 2 ) 


dH( 0, 0, Z) c£_ d z H( 0, 0, Z) 


dZ 


dZ 1 


where Z/(0,0,Z)is given by (6.51). 

Historically, the convection Mach number M c has typically been set equal to 
0-68 Uj/c. (where Uj denotes the mean jet exit velocity) — which would certainly 

be reasonable upstream of the end of the potential core, but would not make 
sense downstream of this point. We therefore generalize this result by setting 


M c = 0.68U c/l / c x (8.3) 

where U c/l denotes the centerline velocity. 

Equation (7.12) can now be summed over the noise producing region of the 
jet to predict the radiated sound. But, as with any acoustic analogy approach, this 
requires information about the mean flow and turbulence statistics, which 
ultimately has to be obtained from measurements of the flow, which limit the 
applicability of the result to a limited range of operating conditions. But it is 
necessary to extrapolate the experimental results beyond the available data base 
in order to extend this approach to a broad range of operating conditions. We do 
this using the experimental results to relate the source model parameters to the 

k-e length, time and velocity scales , k / , y and k V2 .respectively. We then 

use the NASA Glenn version of the NPARC Wind code (Nelson and Power, 

2001 ) to determine the mean flow along with the spatial distribution of the length, 
time and velocity scales within the jet along with the dimensionless time scale 
parameter kU'/e , where the prime is meant to denote differentiation with respect 
to rjn the present context. More specifically, we assume that the length, time 

and velocity scales l x ,l T y x and a 00 that appear in equation (6.48) (directly and 
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implicitly through equations(6.37) and(6.43)) are proportional to the k-e length, 
time and kinetic energy scales * 3/ ^/ ,*/ and (pk) 2 respectively , i.e. we put 

h* c i k3/2 X’ *-**»/. ^ _1 ~ c t%’ % = Ri j kl (y, 0, 0) « C ijkl (pkf (8.4) 

where the dimensionless parameters cj,cy,c T and C ljkl can, in general, depend 
on the dimensionless time scale parameter k\U'\/e. The first three of these are 

determined with reference to the integral length and time scales measured by 
(Bridges & Pobody,1999), and the last is inferred from their mean square 
fluctuating velocity measurements. As noted in section 6, f 1212 depends on the 

autocovariance of the cross velocity components, while the remaining spectral 
components in the acoustic formulas depend only on the covariances of the 
squares of velocity components. We, therefore, expect the parameters in the 
former to scale somewhat differently than those in the latter. 

Since RANS-type models are based on the assumption that the turbulence is 
close to equilibrium, we expect the dimensionless parameters cj,cy,c T and C ljkl 

to be relatively constant in the initial mixing layers and fully developed region 
further downsteam where the flow evolves fairly slowly in the streamwise 
direction (Pope, 2000, p.362 & p. 365) — but not necessarily in the transition 
region, where the flow can be fairly far from equilibrium (see figure 2). We 
account for variations in the latter region by allowing these parameters to depend 
on£|i/'|/e. A comparison of the experimentally measured length, time and kinetic 

energy scales with the corresponding RANS-derived quantities showed that the 
data could be reasonably well represented by using constant values of 
C h C T ,C ijkl everywhere in the flow, but that C T would have to depend onk\U'\/e 

in the transition region. Figure 5 shows results obtained by plotting the ratio of 
the experimentally measured integral time scale obtained from auto-correlations 
of the streamwise velocity fluctuations to the k/s RANS based scale vs. k\U'\/e 

for a Mj = 0.9 jet. These scales are directly proportional to one another in the 
fully developed region and in the initial shear layer, which maps into a single 
point in this plot because k\U'\/e turns out to be nearly constant there. The time 

scale A. 1 can, therefore, be made to transition from its initial shear layer behavior 
to its behavior in the fully developed region by putting C T =C T (0) F(k\U'\/e) . We 

use C T (0) = 0.7 in v f / 1212 C/ 0) = 0.3 1 for all the other spectral components. 

As noted in section 6.2, the individual a 00 values are equal to the density- 

weighted averages of the appropriate velocity components, and it turns out that 
the ones that actually appear in the final formulas only involve squares of these 
components. Lacking any information about the 4 th order correlations we assume 
quasi-normality and infer from equation (6.17) that the dimensionless c ijkl can be 

express expressed entirely in terms of these quantities, or more specifically that 
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C =2 
Mill ^ 


1 2 


? C2222 C 22 „ 2 1 v : 


• c = c = c 

9 M 122 Ml33 M212 


= 2 Vj 2 v' 2 , where 


v' = {vpVjjVj} . The last relation follows from equation(6.19) which implies that 

[Rn (y,0, 0)] 2 =Rn (y,0,0)xR% (j, 0,0) since R xlu (y,0, 0)and R™ 22 (y,0,0) must 

be equal by definition. The mean square velocity components are now assumed 
to be proportional to the turbulence kinetic energy distribution with the best 

overall fit to the data corresponding to the relative proportionality v( 2 = 0.8k and 


v' 2 2 = v' 2 = 0.6k . The anisotropy in the turbulent kinetic energy has a strong effect 
on the 90° degree spectral shape because there are strong cancellations among 
the various components of the spectra x Y*. kl of the generalized Reynolds stress 


autocovariance tensor that contribute to the spectra cD* w appearing in (6.27). 


The (Bridges/Pobody ,1999) data suggest that C x = 0.96, C T = 0.144, for 
v F 1212 and C, =0.58,C r =0.186 for all other spectral components would be 
reasonable values of the remaining parameters in (8.4). It also seems consistent 
to put a = 0.2 , a l 0 = 0.972 a 0 0 for T^and a = 0.34 , a, 0 = 0.54 for all other 

components of (8.2). The remaining unknowns in (7.12), including quantities 
related to the critical layer, depend on the mean flow and can, therefore be 
determined directly from the RANS calculation. 

The WIND code was used to obtain RANS solutions for cold jets with 
acoustic Mach numbers Mj = Ujlc x of 0.50, 0.90, and 1.4--which span the 

range of current interest--and upstream nozzle conditions specified in terms of 
appropriate plenum temperature and pressure ratios. Equations (7.12) to (7.17) 
were then used along with the parameter choices described above to calculate 
their far-field acoustic spectra on the arc x/Dj= 100. 

Figure 6 is a comparison of the results with narrow band data obtained from 
the NASA Glenn SHJAR rig (Khavaran et. al 2002) for the same upstream 
conditions. No vertical adjustments were made to obtain the best fit. Atmospheric 
attenuation was removed from the measurements in order to make lossless 
comparisons with the predictions. The overall agreement appears to be quite 
good but there is a tendency to under predict the higher frequency components 
of the 90° spectrum in the supersonic case. This is because the flow is not 
correctly expanded in this case and the present analysis does not account for the 
resulting shock associated noise (shown cross hatched in Figure 6c). Another 
significant discrepancy is in the prediction of the 0 = 30° spectrum for M d = 0.9 jet 

(figure 6e). This may, in part, be due to the breakdown of the weakly non-parallel 
flow approximation in the region near the end of the potential core - which is an 
important source region for this transonic Mach number. Notice that the spectra 
tend to be narrower and more highly peaked at 0 = 30° than aid = 90° , especially 
at the higher Mach numbers. Figure 7 is a comparison of the present result with 
the best previous results obtained from the Glenn JeNo code, which is based on 
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a parallel model (with an ad hock assumption used to eliminate the critical layer 
singularity) and does not account for the strong retarded time variations 
produced by the long range correlations. 

9. Discussion 

The acoustic predictions of the last section were based on the special case 
(7.12) of the more general result(6.27) which, as pointed out at the end of section 
6.1, is composed of three distinct terms. The first of these (referred to as term / ) 
is the product involving the first term in square brackets. It is not hard to see that 
it exhibits the usual quadrupole- like behavior at low frequencies, i.e. it behaves 
like co 4 as<z> — > 0. The second term (referred to as term// ) is the product 
involving <f>* 212 . It exhibits dipole- like behavior at low frequencies, i.e. it behaves 

like or as® — > 0. It should not, however, be interpreted as an actual dipole 
source, since it is now known (Goldstein (1975), Balsa (1977), Musafir (1992) 
and Afsar et al (2007)) that the mean flow can raise the efficiency of certain 
components of the quadrupole source from co 4 to co 2 as co -» 0. Afsar et al (2007) 
even show that, these components exhibit the same inverse six Doppler factor 
singularity as term//, which as noted above , corresponds to the strongest 
critical layer singularity in the present result . So, while the Goldstein (1975) and 
Afsar et al (2007) result shows that term II is the dominant small angle noise 
source at low frequencies, the present result shows that this term is the dominant 
small angle supersonic noise source at all frequencies. The remaining term does 
not appear to make a significant contribution to the sound field for any of the 
conditions considered and will, therefore, not be discussed any further. 

Term //only contributes to the 90° sound field at extremely low frequencies, but 
becomes increasingly important at small angles as the acoustic Mach number 
increases. Each of the three terms acts like a statistically independent sound 
source. Their individual contributions to the overall 30° spectrum are shown in 
figure 8. Term III is too small to show up on the figure. Notice that Term II is 
only dominant at relatively low frequencies when M , is equal to 0.5, but becomes 

dominant at all frequencies when the Mach number is supersonic. This occurs 
because the six powers of the Doppler factor that appears in the first member of 
this term causes it to be very large. 

Figure 9 shows the contribution of various axial slices of the jet to the 
overall spectra for the supersonic case. Notice that these point-wise spectra 
tend to be much narrower at 9 = 30° than at 0 = 90° -especially in the upstream 
region. This is, in part, due to the fact that much of the small angle sound is 
generated in the critical layer while none of the 90° sound is generated there. 

For acoustic Mach number is close to unity, much of the small angle radiation 
is generated near the end of the potential core -which may explain the tendency 
to under predict the peak frequency of the 30° spectrum in the M s =0.9 case, 

since the weakly non-parallel flow assumption tends to breakdown in this region. 
But this may also be caused by the break down in the quasi-equilibrium 
assumption (on which the RANS modeling is based) in this region. These 
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difficulties can probably be overcome by using a hybrid RANS/LES approach of 
the type described by Goldstein (2006), but that would greatly increase the 
expense of the computation. 

The spectra of each of the three components of (6.27) is determined 
from(8.2), but with different coefficients for each component. The 2nd term in the 
curly brackets of this latter equation is associated with source non-compactness 
effects and does not contribute to the 90° sound field. It is also relatively 
insignificant at all angles when M becomes small but becomes the dominant 
source of the small angle sound atM y = 1.4 . This term introduces an or factor 

which causes term //to behave more like term /at high Mach numbers. 

Computations were only carried out for subsonic convection Mach numbers 
which correspond to the relatively low (<1 .5) acoustic Mach numbers where non- 
linear propagation effects and Mach wave radiation are believed to be 
unimportant. The instability wave contribution to the Green’s function is then 
negligibly small and the convective amplification of the sound sources turns out 
to be relatively insignificant compared to the very strong critical layer effects 
produced by the mean flow Doppler factor. Fortunately, this also corresponds to 
the Mach number range of most technological interest. 

The critical layer only appears when the mean flow acoustic Mach number is 
supersonic and the observation angle (as measured from the downstream axis) 
is fairly small (<45°). It gradually moves inboard from the nozzle lip line with 
increasing downstream distance until it reaches a point beyond the end of the 
potential core where it quickly moves onto the jet axis and suddenly disappears. 
But the propagator is still very close to being singular, and can consequently be 
relatively large, for a significant distance downstream of this point. It was, 
therefore, necessary to construct a new inner solution for this region as well — 
especially since much of the small angle sound field is actually generated in this 
relatively localized portion of the jet when the acoustic Mach number is close to 
unity. This tends to decrease the width of the corresponding acoustic spectrum 
because the band width of the turbulent source tends to be relatively small in 
each streamwise slice of the jet as compared to the relatively large frequency 
variation that occurs over the entire jet. 

The source distribution for the larger angle acoustic radiation tends to be 
more evenly distributed over the jet (since no singular layer occurs at these 
angles). This contributes to making the large angle acoustic spectrum more 
broadband than the small angle spectrum when the jet acoustic Mach number is 
close to one, but the relative spectral width of the individual slice contributions is 
probably the more important factor here. 

The upstream critical layers move into the high turbulence region along the 
nozzle lip with increasing jet Mach number. This causes more of small angle 
radiation will be generated in these layers, which generate narrower spectra than 
the reminder of the jet when the correct inner solutions are used in these 
regions. This partially explains why the width of the small angle spectra 
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decreases with jet Mach number. It is, of course, well known that Doppler 
frequency shift and mean-flow refraction effects also play an important role here. 
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Appendix A Mean Flow 

As in Goldstein & Leib (2005), we suppose that all lengths have been normalized 
by some characteristic cross flow dimension of the jet (say its diameter, in the 
case of a round jet) and all velocities by some appropriate characteristic 
streamwise velocity with similar obvious normalization for the density, pressure 
and temperature. Then the mean flow variables and stresses should expand like 


U(Y,y T ) + £U (l) (Y,y T ) + ... 

(A. 1 ) 

eV(Y,y T ) + s 2 V (l) (Y,y T ) + ... 

(A.2) 


p e =P(Y,y r ) + eP {l) (Y,y T ) + s 2 P (2) (Y,y T ) + ... (A.3) 


p = R(Y,y T ) + sR (x \Y,y T ) + ... (A.4) 

and 




= -pv t 



+ (j , = t^) + st}:> + + 


to 


■^( 2 ) 


vj 


VJ 


y “J 




(A.5) 


where 

Y = s y x (A. 6) 

denotes the slow streamwise “source” variable, and 


y T ={t 2 ,t 3 } 


(A.7) 
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V T = {v 2 ,v 3 } 


(A.8) 


denote cross flow variables. 

Substituting these into the mean flow equations (2.7) to(2.9) (written in terms 
of the “source” variable y ) and assuming that the generalized Reynolds stresses 
vanish in the free stream, shows that the result will only balance if 

Tip = Tpp = Tip = 0 for y=2,3 (A.9) 

( P _ T jP ) = TT~ (nosumony *1 = 2,3) (A. 10) 

oy } v ’ ay i 

and 


^_( p (1) - 7 ^ 1) ) = (no sum ony ^ / = 2,3) . (A.11) 


The lowest order mean flow equations become 


DR= 0 


D RU = S 


dP (2) 8 


(i) , 01 J‘ 
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T’i 1 ) 

1 4j 
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+ UTP +V,Tp >) +- 


Typ’ + t t) 


for/,/ = 2,3, 


(A. 12) 
(A. 13) 

(A. 14) 


(A. 15) 


where the operator £> 0 is now given by 


and we have put 


ps ps 

D 0 = —U + ——V, for j = 2,3 

0 BY By j 1 


(A. 16) 
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S=^(t!?-P) + 1 ?-TI;\ fory = 2, 3 


Equations(2.21 ), (A.3) and (A. 5) imply that 


s, =s,( p+epl ")-( T !°' +eT P) +0 ( s2 ) 


and it therefore follows from (A. 10) and (A.1 1 ) that 


= -^,s+0( £ ! ). 


Appendix B Details of Critical Layer Solution 

Equation (4.35) can be written as 


K-uM ^--2a ! y ± u( d -^\ +K u: si 


8Y J[8y 


or since the chain rule implies that 


U.^JdU] dU(y„(Y),Y) p 
c 8Y { 8Y J dY 


— -2ay L jg" 4 =U' c gl 


And it, therefore, follows that 
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(B.4) 


which implies that g 5 a 4 °=0 in (4.27). Inserting (4.41) into (4.33) shows that 
g° 4 satisfies 
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which, in view of (4.44) possesses the solution 
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y± 


that matches with the 2nd equation(4.27). And inserting this along with 
equations(4.43) and (4.45) into equation(4.34) and using (4.29) and (4.39) shows 
that 
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which, in view of the 2 nd equation(4.27), implies that 


(B.7) 
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Appendix C Uniformly Valid Formulas for Propagator 
Components 

In the “outer region” where the source point, y T , is outside of the critical layer, 
equations(4.21 ),(5.20) and (5.32) imply that 
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v C 0O 


0,®| y T 


(C.1) 
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(C.2) 


Gu = ~ 


i cos 9 


2 cm 


d & 0 + d 




(1-Mcos0) 2 <tyi dy t (1-Mcos0) 2 


+ 1 


,y-\dU 


2co dy t (l-Mcos9) 


co 


(1 -Mcos9) 


cos 




d k/Xc 


(y-1) dM~ 



c x dy f (1-M cos9) 2 * dy t 


i = 2,3 (C.3) 


where the Doppler weighted cosine function C is defined by 


2 c 2 cos 7 6 


C 1 = 


cl [M(j r )cos0-l]' 


(C.4) 


To obtain the corresponding formulas in the “inner region”, notice that 
equations (4.54) and (5.22) imply 
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a[y,k,a\y T )= (/ 0 [y,k,a\y T )\ (C.6) 

y j_ 

and use the coordinate transform 

yi,y 2 >y3^>yi>y±> s ( c - 7 ) 

together with equations (4.48) and (4.51 ) to obtain 

o, =T ll =^- = -a.[l + w'(aJ 1 )], (C.8) 


G„ =- 2a ' , ‘' . *✓<«„,, ), 

<=(g;) 


where the prime denotes differentiation with respect to the entire argument and 



It therefore follows that 
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Then since equations (4.48) to (4.50) show that 
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it follows from the coordinate transform (C.7) that 
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to the lowest order of approximation. And similarly that 
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Using the product rule (see discussion at the end of section 4), to combine 
equations (C.1) to (C.17) and retaining only the lowest order approximation within 
the critical layer, we obtain the uniformly valid approximations to the sums and 
products 
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of the dominant G.. components, where z 0 and C c are given by (6.15) and (6.16). 

Appendix D Isotropic Source Model 

Inserting equations(6.1 1 ) and (6.12) into(5.31 ) yields 
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where G, is defined by(5.32). Then since 
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And, in view of (C.18) to(C.20), it follows that 

Gij { x \yT)^ij (• X: |j ; 7 T ) _ (x|j r )| =2 |G 21 | + |G 31 | + 1^23 I “ ^(^ 22^11 + ^33^11 + ^22^33 ) 


= 2 


|g 2 ,| 2 + |g 3i | 2 + |g 23 | 2 +|g „| 2 -»(g„g* +g 22 g 3 ' 3 ) 


(D.4) 


= 2 


|g ,,| 2 +|g 3I | ! +|g 23 | 2 -®g, 2 g ; 3 +®c t (c; -i)[i+w '( Zo )] 2 


a 


9 ,— cos 9, col j r 

V C oo J 


and, therefore that 

7 »(*W 

44 »C c (c;-l)® 


r 2n^ 

27TCO . ^ 



sin 0 < 

(T) ^ 

1 mi 

v X ) 

Ao i 

V 


y;— cos 0 , — V y . r S, co(l-M c cosG) 

'no ^ 


A 


f 


1212 


y; — cos 0, — V S, to (l - M c cos 0 ) 

V ^00 ^00 


A 




cp, — cos 0 ,co|j r [l + w'(z 0 )] 


V C oo 


+4 

(D.5) 


G 2 1 + Cr 31 + (|(j r J — /^eG 11 G 


'23 


22 w 33 


Z' 


o : 


1212 


y;— cos0, —V yT S, co(l -M c cos0) 

V ^ 00 Ao 


A 


where C is given by(6.16). And, it now follows from (C.21) and(C.22) that 
equation (6.14) holds. 

Appendix E Axisymmetric source Model 

Goldstein & Rosenbaum (1973 a) show that 

H i] G kl‘ikj l = H » G lAn ~(n n Gl + H„G’ 22 )(I, m -I^ + H'jGyQy (E.1) 

with I jkl defined by equation (6.23) Q n = I nn -I 1212 , Q n = Q n = I nn +/ 1122 , 
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022 = 033 = 0 222 ~ him 023 = 0 222 ~ 0323 for any 2 nd order symmetric tensors G t] 
and H tj and any axisymmetric tensorial function R™ (t) of %. It, therefore, 
follows from equations (D.2) and(D.3) that 
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It, therefore follows from (C.18) to (C.22) and (6.16) that equation (6.24)holds. 
Taking H tj = 5 , y shows that 
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And taking H tj = G* = 6 y shows that 
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Appendix F Details of Spectral Component Model 

Integrating equation(6.34) by parts shows that 
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where D k , D m , and D k/ are defined by equations(6.32). Then since (Campbell & 
Foster, 1942, pill, #867) 
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e mT ~ x dr = (F.2) 
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where ^is a modified Bessel function of the second kind. (Campbell & Foster, 

1942, pi 25, #917.8) and equation (6.30) show that equation (F.1) can be written 
as 
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which can now be rewritten as equation (6.35) . 

Appendix G Axisymmetric Jet 

When expressed in polar coordinates, the coefficients in equations(6.14) 
and ( 6 . 27 ) become 
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where C is defined by (C.4). Since 
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equation (G.2) becomes 
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Integrating these over the circumferential direction and using equation(7.1 1 ) 
shows that 
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Inserting these into equations (4.21 )and (6.27) shows that (7.12) holds. 
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Figure 1 Critical layer surface 
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Figure 2 Flow Configuration 


Outer Region 



Figure 3 Critical Layer Structure in Supersonic Jet 
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Figure 4 Comparison of the functional form (6.32) with experimental data. 












Part e)Mj =0.9, 0 = 30° 



Part f) Mj = 1 .4 0 = 30° 

Figure 6 Comparison of Predictions with data 


7i 





Figure 7. Comparison with best previous results obtained from Glenn JeNo code 
(Khavaran, Bridges,. & Georgiadis, N., 2005) 
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Part (c) Mj =0.5 

Figure 8. Contribution of individual terms to overall spectra at 0 = 30° 


0.0 < =X/D < = 2.0 
2.0 < =X/D < = 4.0 



Part a) 9 = 30° 
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0.0 < =X/D < = 2.0 
2.0 < =X/D < = 4.0 



Part b) 9 = 90° 

Figure 9 Contribution of individual slices to overall spectrum for Mj =1.4 jet 
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